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SUMMARY 


V 

*  new  two  parameter  family  of  Ilf,  length  distribution, 
U  Pr'""Ced  uhlch  “  derived  from  .  model  for  fatigue,  Thi. 
derivation  follow,  from  consideration,  of  renewal  theory  for 
the  number  of  cycle,  needed  to  force  a  fstlgue  creek  exten.ion 
to  exceed  a  critical  value.  Some  clo.ure  propertle.  of  thi. 
family  are  given  and  some  comparison,  made  with  other  families 
such  aa  the  lognormal  which  have  been  prevloualy  used  In 
fatigue  studies. 
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2.  A  MODEL  FOR  THE  DISTRIBUTION  OF  LIFE 

We  propose  a  basic  framework  and  notation  which  is  similar  to  that 

used  previously  in  [  3  ] .  We  consider  only  standardized  material  specimens 

which  are  subjected  to  fluctuating  stresses  by  a  periodic  loading.  By  a 

load  (or  load  oscillation)  we  mean  a  continuous  unimodal  function  on  the 

unit  Interval,  the  value  of  which  at  any  time  gives  the  stress  imposed 

by  the  deflection  of  the  specimen.  Let  i^,^****  be  the  sequence  of 

th 

loads  which  are  to  be  applied  at  each  oscillation  so  that  at  the  1 
oscillation  load  9.^  is  imposed.  We  suppose  that  the  loading  is  cyclic 
in  the  sense  that  for  some  m  >  1  and  all  i"l,...,m 

£jm+i  "  *knri-i  for  a11  ^  **  k  (2.1.1) 

and  the  loading  is  continuous  so  that  for  all  i»l,2,... 

£i+l(0)  "  £i(1)'  (2.1.2) 

8 1 

Hence  the  (j+1)  cycle  is  the  loading  * 

We  assume  that  fatigue  failure  is  due  to  the  initiation,  growth  and 
ultimate  extension  of  a  dominant  crack.  At  each  oscillation  this  crack  is 
extended  by  some  amount  which  is  a  random  functiun  due  to  the  variation  in 
the  material,  the  magnitude  of  the  imposed  stress  and  a  certain  number  of 
the  prior  loads  and  perhaps  the  actual  crack  extensions  caused  by  the 
prior  loads  in  that  cycle. 

Thus  we  now  make  our  first  assumption 

1°  The  incremental  crack  extension  following  the  application 
of  the  1th  oscillation  is  a  random  variable  with  a  distribution  which 
depends  upon  all  and  only  the  loads  and  actual  crack  extensions  which 
have  preceded  it  in  that  cycle. 
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1.  INTRODUCTION 


It  is  well  known  that  for  the  amount  of  fatigue  data  which  can 
usually  be  obtained  almost  any  two  dimensional  parametric  family  of 
distributions  can  be  made  to  fit  reasonably  well.  In  fact,  in  the 
region  of  central  tendency  the  lognormal,  the  Weibull,  the  Gamma, 
etc.,  can  all  be  fitted  by  parametric  estimation  and  because  of  the 
relatively  small  sample  sizes  hardly  any  can  be  rejected  by,  say  a 
Chi-square  Goodness  of  Fit  test.  However,  when  it  becomes  a  question 
of  predicting  the  "safe  life"  say  the  one  thousandth  percentile, 
there  is  a  wide  discrepancy  between  these  models. 

For  this  reason  a  family  of  distributions  which  is  obtained  from 
considerations  of  the  basic  characteristics  of  the  fatigue  process  should 
be  more  persuasive  in  its  implications  than  any  ad  hoc  family  chosen 
for  extraneous  reasons.  In  this  paper  we  derive,  using  some  elementary 
renewal  theory,  a  two  parameter  family  of  nonnegative  random  variables 
as  an  idealization  of  the  number  of  cycles  necessary  to  force  a  fatigue 
crack  to  grow  to  a  critical  value.  We  then  examine  some  of  its 
relevant  properties. 
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This  assumption  can  be  plausibly  held  for  ground-air-ground 
cycles  in  aeronautical  fatigue  studies  and  other  such  applications. 

s  t 

The  crack  extension  during  the  (j+1)  cycle  is 


j+1  jnH-1 


+  •••+  X 


jtrrhn 


for 


j-0,1 


t 


where  X 
the  load 


jm+i 

l. 


is  the  (possibly  microscopic)  crack  extension  following 
applied  in  the  i*"*1  oscillation  of  the  (j+l)St  cycle. 


It  follows  from  Assumption  1°,  regardless  of  how  much  dependence 
exists  between  the  successive  random  extensions  per  oscillation  in  each 
cycle,  that  the  random  total  crack  extensions  per  cycle  are  independent. 

Thus  we  could  formally  make  a  second  assumption 

2°  The  total  crack  extension  due  to  the  jth  cycle  is  a 

random  variable  with  mean  p  and  variance  for  all  j«l,2,...  . 


Our  notation  will  be 


W 

n 


n 


which  has  distribution  function 


H  (w)  ■  P[W  <  wj,  for  n-1,2 . 

n  n 

It  follows  from  elementary  probability  considerations,  see  p.  189, 
[6],  that  the  distribution  of  C,  the  number  of  such  cycles  until 
failure,  in  the  case  failure  is  defined  as  the  crack  length  exceeding 
some  fixed  critical  length  m  for  the  first  time,  is 
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TOKWV!  a 
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P[C  £  n]  -  1  -  H  (oj). 
_  n 


(2.2) 


Now 


P [C  £  n]  -  P 


Is- 

Li-lo 


i -v 


>  w~ny 

}/n  o/n 


(2.3) 


and  since  the  are  independent  and  identically  distributed  random 

variables,  using  the  central  limit  theorem  for  n  large,  Equation  (2.3) 
can  be  approximated  by  the  normal  distribution.  Let  91  be  the 
distribution  function  of  the  standard  normal  variate  with  zero  mean 
and  unit  variance  defined  for  -®  <  y  <  ®  by 


ry  1 

9Ky)  -  /  -zz 

v2v 


e-t2^2dt. 


(2.4) 


We  now  make  explicit  this  assumption  of  equality 

3°  the  distribution  of  C  from  (2.2)  is 

P[C  <.„].!  -  9/4^  -  • 

\/n  0/  \  /no/ 

We  are  cognizant  of  the  fact  that  Assumption  3°  takes  as  exact, 
for  the  applications  intended,  an  approximate  equality.  Such  an  assumption 
can  only  be  justified  on  physical  grounds,  not  mathematical. 

Firstly,  if  we  suppose  that  we  are  dealing  with  a  long  cycle  of 
oscillations  which  is  as  complex  as  the  ground-air-ground  cycle  in 
aeronautics  it  might  be  reasonable  to  assume  that  each  cycle  itself 
consists  of  a  large  number  of  distinct  phases  of  loading.  Even  though 
the  total  crack  extension  per  cycle  is  the  sum  of  random  variables 
which  are  not  necessarily  either  independent  or  identically  distributed, 
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the  number  of  summands  might  be  sufficiently  large  to  make  it  reasonable 
to  assume  that  the  total  crack  extension  per  cycle,  namely  Yj ,  is  itself 
approximately  normal. 

Secondly,  because  in  the  case  mentioned  above,  fatigue  is  measured 
in  the  thousands  of  cycles  and  in  many  other  applications  in  the  millions 
of  cycles,  it  is  felt  that  such  large  numbers  of  random  variables  make 
3°  an  eminently  practical  assumption  especially  in  view  of  its  simple 
analytic  form. 

This  distributional  form,  parameterized  differently,  has  been 
previously  obtained  by  Freudenthal  and  Shinozuka  in  [7].  In  that  report, 
which  was  unpublished,  only  a  heuristic  derivation  based  on  engineering 
considerations  was  presented  and  then  ad  hoc  fitting  procedures  were 
used  to  substantiate  the  validity  of  that  model  with  several  sets  of 
fatigue  data. 

As  a  possible  alternative  one  could  postulate  that  the  distribution 
of  the  crack  extension  might  be  different  for  the  earlier  cycles  than  It 
would  be  for  the  later  ones.  One  of  the  first  suppositions  might  be  to 
make  the  distribution  of  the  extension  per  cycle  depend  upon  the  size  of 
the  crack  at  the  start  of  the  cycle.  One  such  assumption  and  the  result¬ 
ing  distribution  of  the  total  crack  length  at  the  end  of  n  cycles  is 
presented  in  Section  4,  We  do  not  make  that  assumption  here  but  instead, 
for  reasons  of  simplicity,  proceed  with  the  study  of  the  implications 
of  the  one  above. 


We  write 


» 


a 


6 


w_ 

y 


(2.5) 
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and  replace  n  by  the  nonnegative  real  variable  t  >  0.  If  we  now 
denote  the  continuous  extension  of  the  discrete  random  variable  C  by 
T,  a  continuous  nonnegative  random  variable,  then  it  follows  by  3° 
and  (2.5)  that  T  has  the  life  distribution 


F(t:a,fO  -  £(*78)1  for  t  >  0  (2.6) 

a 

where 


a  >  0,  8  >  0 


and 

£ (t)  -  t*  -  t"*5. 


(2.7) 


This  two-parameter  family  of  distributions  is  a  plausible  model  for 
the  distribution  of  fatigue  life.  The  set  of  all  life  distributions 
of  the  form  (2.6)  for  a,6  >  0  will  be  denoted  by  £7  and  most  of 
this  paper  will  be  devoted  to  the  study  of  its  properties.  We  shall 
also  refer  to  the  law  which  has  the  distribution  (2.6)  with  the 
somewhat  shorter  notation  F(a,8). 

Of  course,  there  are  motivational  derivations  for  other 
distributions  as  well.  The  Weibull  distribution,  which  is  well 
known  for  its  applications  as  a  life  length  distribution  and  for 
fatigue  life  in  particular,  is  obtained  as  a  special  case  of  the 
extreme  value  distributions,  see  p.  302,  [8].  The  Gamma  family 
has  also  been  obtained  as  a  distribution  of  life  by  utilizing  a 
model  of  a  bundle  of  strands  which  are  supporting  a  tensile  load, 
see  [A]. 

It  is  instructive  to  make  a  comparison  between  the  derivation 
of  the  family  of  distributions  &  and  an  appropriate  adaptation  of 
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the  classical  heuristic  argument,  found  for  example  p.  219,  Cramdr  [5], 
as  it  might  be  used  to  obtain  the  lognormal  distribution  of  the  time 
until  failure  in  fatigue.  This  argument  is  presented,  along  with 
derivations  of  other  distributions,  by  Parzen  in  [9].  We  should  also 
mention  that  damage  accumulation  due  to  fatigue  was  first  treated  as  a 
renewal  problem  in  that  paper. 


m 
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3.  SOME  PROPERTIES  OF  9 

Let  T  have  the  distribution  defined  in  (2.6).  Note  that  T 

is  a  two-parameter  random  variable  with  8  as  a  location  parameter 

since  it  is  the  median  of  the  distribution.  (We  show  later  that  8 

is  neither  the  mode  nor  the  mean.)  Notice  also  that  a  is  a  shape 

parameter  and  8  a  scale  parameter.  As  we  have  seen  ^  £(T/8)  is 

a  standard  normal  random  variable  with  mean  zero  and  unit  variance. 

2 

If  we  let  X  be  9J(0,  “)  we  see  that,  in  distribution, 

2X  -  £(T/8).  (3.1) 

If  we  define  the  function  tp  by 

<p(x)  -  £  1(2x)  for  all  real  x,  (3.1.5) 

then 

T  -  8*(X).  (3.2) 


From  elementary  algebra  we  find  that 


*(x)  -  [p(x)]2 


where 


o 


p(x)  -  X  + 


/x2+l 


So  by  (3.2) 


T  -  6[1  +  2X2  +  2X/1+X2  ) 


where  X  is  9?(0 ,  — ) .  Hence  we  have  immediately 

a2 

E(T)  •«  6(1  +  y) 


(3.3) 


(3.4) 


(3.5) 


(3.5.1) 
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E(T2)  -  02(1  +  2a2  +  •—) 

2  5a2 

var(T)  -  (a0r(l  +  ^L-)  (3.5.2) 

and  we  note  that  for  fixed  a  the  variance  of  T  increases  as  the  scale 
parameter  (median)  0  increases.  This  is  not  true  for  the  lognormal 
distribution  but  empirical  evidence  shows  that  it  is  true  for  the  fatigue 
lives  themselves. 

By  noting  that 

-  -x  +  £7.  p(-x)  (3.5.3) 

we  see  that,  whenever  -X  has  the  same  distribution  as  X,  we  have  by 
(3.3),  in  distribution, 

“d  ^o‘*(x)-  (3-6) 

Thus  there  follows  immediately  from  (3.2)  the 

Theorem  3.1.  If  T  has  a  fatigue  life  distribution  F(a,0),  in 
9  then  y  has  a  distribution  in  9  given  by  F(a,  j) .  Moreover,  for  any 
real  a  >  0,  the  random  variable  aT  has  a  distribution  in  9  given  by 
F(a,a0) . 

It  is  known  that  every  random  variable  with  distribution  defined  by 
(2.2)  for  which  the  Yj  are  nonnegative  and  have  densities  which  are  P61ya 
frequency  functions  of  order  2,  has  an  increasing  failure  rate,  see  [1]. 
While  our  random  variable  T  does  not  have  this  property,  its  average 
failure  rate  is  nearly  nondecreasing.  Specifically,  it  has  the  properties 
of  our  subsequent  remarks  (2.2)  and  (2.3).  The  class  of  failure  rates 
which  are  nondecreasing  on  the  average  has  been  studied  in  [2]  and  was 
shown  there  to  have  some  closure  properties  which  are  physically  plausible. 


Will  t  nw  WWVW 


-10- 


Remark  3.2.  T  has  an  average  failure  rate  which  approaches 
a  positive  constant. 

Proof.  Without  loss  of  generality,  take  a  ■  6  ■  1  and  define 
Q(t)  ■  -£n{l  -9UC(t)]},  h(t)-^1. 

Writing  Mill's  ratio 


M(0  -  ui  -  9K01/WU) 

then  applying  L'Hospital's  rule  to  h(t)  for  t  ■*  00  and  using  (A. 1.1) 
of  the  next  section  we  obtain 


h(t) 


21*  IS.lt)]  „ 
l-9JU(t)] 


C'(t) 


ii-t;2) 

2M[£(t)]  • 


We  know  that  as  £  ■+  ®  we  have 


M(£)  =  1  "  -T  +  -T+  0a"5). 

r  r 

Hence  lim  h(t)  ■  j  but  Q(l)  *  in  2  >  y  .  This  proves  our  contention 
that  h  does  not  always  increase  but  approaches  a  positive  constant. || 

Remark  3.3.  Actual  numerical  calculation  shows  that  the  average 
failure  rate  h  decreases  slowly  for  t  >  1.64,  as  illustrated  in  Figure  1. 

Surprisingly  enough  we  regard  this  particular  feature  of  the 
distribution  as  a  virtue  since  it  is  in  agreement  with  the  observed 
facts.  Typical  fatigue  data,  such  as  given  in  [7],  show  that  we  cannot 
assume  that  the  failure  rate  is  always  increasing.  Such  data  can  be 
explained  as  being  the  mixture  of  distributions  each  with  increasing 
failure  rate  such  that  the  combined  failure  rate  decreases  for  a  time. 


*  .  *  -■J7+*+T* r#. 
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However,  almost  any  model  which  is  the  mixture  of  distributions,  even 
two  Weibull  distributions  with  the  same  shape  parameter  but  different 
scale  parameters,  would  necessarily  be  so  cumbersome  mathematically  that 
the  statistical  problems  which  arise  would  be  difficult  analytically. 


Figure  1 


Graphs  of  the  density  f,  the  distribution  F,  the  hazard  rate 
and  the  average  hazard  rate  h  where  for  t  >  0 

f(t>  -  ~=T"  exp{-  f  -  7T+  1} 

2t/2it 

F(t)  -  J  f (x)dx,  q(t)  -  » 


a.  k. 


h(t) 


u 


q(x)dx. 


*»Hgl 
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4.  A  SECOND  MODEL  FOR  THE  DISTRIBUTION  OF  LIFE 

In  this  section  we  replace  the  Assumption  2*  by  the  alternative 

one 

6 1 

2'  The  crack  extension.  during  the  (n+1)  cycle,  given 

that  the  total  crack  length  was  s  at  the  start  of  the  cycle, 

is  a  normal  random  variable  with  mean  u  +  6s,  for  some  constant 

2 

6^0,  and  variance  o  for  each  nal,2,...  . 

By  formula  (2.2),  to  find  the  distribution  of  C  it  is  necessary  and 

sufficient  that  we  find  the  distribution  of  W  ,  called  H  . 

n  n 

By  assumption 

HVjiy  I  "n“sl  ' 

and  for  n-1 

H^y)  -»  <  y  <  co.  (4.1) 

Now  by  definition,  setting  y  -  x  -  s 

W*>  ■/  (4'2> 


We  can  now  prove  the 

Theorem  4. 1 .  The  total  crack  length  W^,  at  the  end  of  the 

th  2 

n  cycle,  is  normal  with  mean  u  and  variance  o  where 


u  -  *  [<l+6)n-l], 
n  6 


2  2  (1+6) 2n-l 

o-o  - - -  . 

(1+6) 2-l 


(4.3) 


Proof  by  induction.  The  statement  is  true  for  n*l.  Assume  it  true 


for  n,  then  from  (4.2) 
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w*>  ■  /»«V>d  s91^ 

'■  n 


where 


u 


1+5 


v 


0 

1+5  * 


But  this  we  recognize  as  the  usual  convolution  of  two  normal  random 
variables  and  hence  the  result  Is 


Wx> 


Hence  by  simplification  we  find  that 

Vl  ‘  "  +  <1+«  V  °n+l  '  °2  +  (l+<,2°n 


and  one  checks  that  formulas  given  do  satisfy  the  recursion  relations 

(4. 3). 1 1 

Strictly  speaking  a  proper  distribution,  analogous  to  that 
obtained  in  (2.6),  which  would  be  the  continuous  extension  of  the 
present  case,  cannot  be  generated  since  P[C  <  00 ]  <  1. 

For,  one  can  see  that 

n  -to 

P(C  <  n]  -W(— — ) 

—  a 

n 

and  as  n  -*  00  we  have  p  /o  -*•  J (l+5)^-l  < 

n  n  5o 

Of  course  this  would  not  be  of  practical  significance,  since  we  must 
restrict  ourselves  to  situations  where  ^-  >  3,  otherwise  2'  would  be 

u  _ 

physically  unrealistic.  Hence  lim  —  >  3/2/5  for  0  <  6  1,  and 

p  n-K»  n 

lim  —  >  3  for  any  6. 
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5.  CONCLUSION 

A  derivation,  based  on  plausible  physical  considerations,  for 
a  family  of  distributions  is,  by  itself,  not  a  conclusive  argument 
that  such  a  particular  family  should  be  used  in  life  studies.  No 
family,  however  reasonable  its  derivation,  can  be  accepted  for  use  in 
fatigue  life  studies  until  it  is  confronted  with  actual  fatigue  data 
obtained  under  various  conditions  and  the  distribution  is  shown  to 
represent  adequately  the  life  lengths  which  are  obtained. 

In  order  to  do  this  one  must  have  the  theory  of  estimation  for 
this  family  completed.  The  derivation  of  parametric  estimators  and  the 
ancillary  computing  formulas  for  this  family  will  be  presented  in  a 
latter  study.  Also  further  studies  of  the  application  of  this 
distribution  to  the  calculation  of  "safe  life"  will  be  made.  Thus  the 
confrontations  of  this  family  with  actual  data  will  be  carried  out,  to 
provide  the  justification  for  this  presentation. 


• ,  «?*•  +&& 
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